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In alternative theories of gravity, designed to produce cosmic acceleration at the current epoch, 
the growth of large scale structure can be modified. We study the potential of upcoming and future 
tomographic surveys such as DES and LSST, with the aid of CMB and supernovae data, to detect 
departures from the growth of cosmic structure expected within General Relativity. We employ 
parametric forms to quantify the potential time- and scale-dependent variation of the effective 
■ gravitational constant, and the differences between the two Newtonian potentials. We then apply 

' the Fisher matrix technique to forecast the errors on the modified growth parameters from galaxy 

clustering, weak lensing, CMB, and their cross-correlations across multiple photometric redshift 
bins. We find that even with conservative assumptions about the data, DES will produce non- 
trivial constraints on modified growth, and that LSST will do significantly better. 
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I. INTRODUCTION 



Ph . . 

I ' Observations strongly favor a universe that has reeently entered a phase of accelerated expansion [1, 2]. This poses a 

2 [ puzzle for modern cosmology as standard General Relativity (GR), applied to a universe which contains only radiation 

. and dust, has difficulties fitting the data. One can view this as evidence for the existence of Dark Energy (DE) - a yet 

' unknown component with a negative equation of state, such as a eosmological constant, A. An alternative explanation 

I— I, could involve modifying GR in a manner that leads to accelerating solutions. Popular examples include the so-called 
f{R) class of models [3-8], Chameleon type scalar-tensor theories [9], the Dvali-Gabadadze-Porrati (DGP) model [10], 

^ [ and models motivated by DGP, such as the recently introduced Degravitation scenario [11]. 

• Although the ACDM model, consisting of GR with A and Cold Dark Matter (CDM), is currently the best fit to 
] the data, it faces some challenges on the theoretical side, such as the coincidence and the fine-tuning problems. On 

— t the other side, alternative theories of gravity have to satisfy the multitude of existing experimental tests passed by 
CO ' GR [12, 13]. This typically requires a degree of fine-tuning which, at best, does not improve on that involved in 
, setting A to the value required by current data. Moreover, once these modified theories are tuned to avoid conflicts 
■ with existing constraints (when it is possible), their predictions for the expansion history of the universe are often 
, identical to that of the ACDM model^[ll, 14-18]. However, this degeneracy is typically broken at the level of 

• eosmological structure formation; indeed, models of modified gravity that closely mimic the eosmological constant 
at the background level can still give significantly different predictions of the growth of structure. The large scale 
structure of the universe therefore offers a promising testing ground for GR and it is important to explore to what 
extent one can detect departures from GR in the growth of structure with present and upcoming eosmological data. 

^ . By definition, the term 'modified gravity' implies that the form of the Einstein-Hilbert action is different from that 
■ - ' of GR, and, as a consequence, the Einstein equations are changed. At the background level, the modifications allow 
for a late-time acceleration, which is typically degenerate with ACDM after the required tuning. However, since the 
equations describing the evolution of eosmological perturbations are modified as well, models with the same expansion 
history as in ACDM can lead to different dynamics for the growth of cosmic structure. 

Structure formation has been studied for f(R) models in [16, 19-21], for Chameleon models in [22], and for the 
DGP model in [23-26]. One way to test the consistency of the ACDM model is to compare values of the eosmological 
parameters extracted from distance measures, such as supernovae magnitudes and baryon acoustic oscillations, to the 
values found from growth measures, such as galaxy counts and weak lensing [27]. One can also introduce general 
parametrizations of the modified evolution of gravitational potentials and matter perturbations [28-32] for the purpose 
of detecting/constraining departures from GR. 



^ The modifications of gravity discussed in this work do not attempt to replace dark matter. We assume existence of CDM. 
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Scalar metric perturbations in the Newtonian gauge are described by two potentials, ^{x,t) and t), which 
correspond to perturbations in the time-time and space-space components of the metric tensor, respectively. In the 
ACDM model, these two potentials are equal during the epoch of structure formation, and their time dependence is set 
by the same scale- independent linear growth function g{a) that describes the growth of matter density perturbations. 
This, generally, is no longer true in theories of modified gravity, where one can have scale-dependent growth patterns. 
The two Newtonian potentials need not be the same, and their dependence on matter perturbations can be different. 
Working in Fourier space, we parametrize the ratio between the Newtonian potentials and the dependence of ^ on 
the matter density perturbation with two time- and scale-dependent functions j(a,k) and fi(a,k). We then study 
the potential of upcoming and future tomographic surveys to constrain departures of these functions from their GR 
values. Dealing with unknown functions implies working with an infinite number of degrees of freedom. Similarly to 
the more widely studied problem of determining the dark energy equation of state, w{z), one can proceed in several 
ways: 

1. One can assume a functional form motivated by a certain class of theories, with a few parameters, and forecast 
the constraints on the parameters. This would tell us the extent to which one can constrain these theories and 
also allow us to reconstruct the shape of the functions based on the chosen form. The results, of course, would 
depend on the choice of the paramctrization. However, they would be good indicators of the power of current 
and upcoming surveys to constrain departures from GR. Moreover, as we will discuss shortly, it is possible to 
employ a paramctrization that accurately represents a broad class of modified theories. This is the approach we 
take in this paper. 

2. Another, non-parametric approach, consists of performing a Principal Component Analysis (PGA) to determine 
the eigenmodes of /i and 7 that can be constrained by data [33, 34]. This method allows one to compare 
different experiments and their combinations according to the relative gain in information about the functions. 
PGA can also point to the "sweet spots" in redshift and scale where data is most sensitive to variations in 
jj, and 7, which can be a useful guide for designing future observing strategies. The PGA method does not 
allow one to reconstruct the shape of the functions from data. However, one can still reproduce the errors 
on parameters of any paramctrization from the eigenvectors and eigenvalues found using PGA [34]. Hence, in 
terms of forecasting the errors, the PGA method can do everything that the first method can, plus the benefits 
mentioned above. It is, however, more demanding computationally, and one needs a criterion for deciding which 
modes are well-constrained. Wc consider the PGA method in a separate publication [35]. 

3. A third approach, which can work for certain estimators of 7, is a direct reconstruction from data. In [28, 36], 
it was proposed to consider the ratio of the peculiar velocity-galaxy correlation with the weak lensing - galaxy 
correlation. In such a ratio, the dependence on the galaxy bias cancels out. Then, since the peculiar velocities 
are determined by the potential 4", while the weak lensing is controlled by $ -I- such ratio, if appropriately 
constructed, would directly probe any difference between $ and 'J. This is a more direct and model-independent 
way of testing GR with the growth of structure than the first two methods. Its power, however, will depend on 
how well future experiments will be able to measure peculiar velocities. Also, while it may allow the extraction 
of 7, it does not directly probe fi. Still, this is a novel and promising method that should be pursued in parallel 
with the first two. 

Our parametric forms for 7(0, A:) and ^{a,k) are analogous to those introduced in [32], and contain a total of five 
parameters. These forms are highly accurate in describing the linear growth in a wide class of scalar-tensor theories 
and also fiexible enough to capture features of modified dynamics in other theories. We use the Fisher matrix technique 
to forecast the errors on the modified growth parameters, along with the standard set of cosmological parameters. We 
represent the data in terms of a set of all possible two-point correlation functions (both auto- and cross-correlations) 
between the galaxy counts, weak lensing shear, and Cosmic Microwave Background (GMB) temperature anisotropy, 
across multiple redshift bins, in addition to the GMB E-mode polarization autocorrelation, and the GMB E-mode 
and temperature cross-correlation. We find that even with a conservative treatment of data, such as using only the 
modes that are well within the linear regime, the upcoming and future tomographic surveys, like the Dark Energy 
Survey (DES) and Large Synoptic Survey Telescope (LSST), in combination with GMB and future SuperNovae (SNc) 
luminosity-distance data, will be able to produce non-trivial constraints on all five parameters. We also show that 
LSST will significantly improve our ability to test GR, compared to DES. 

The rest of the paper is organized as follows. In Sec. II, we motivate our choice of parametric forms for 7(0, k) 
and n{a,k). In Sec. HI, we describe the set of observables expected from tomographic weak lensing surveys, their 
dependence on the underlying gravitational potentials, and the Fisher matrix technique we used to forecast the 
parameter errors. We also describe the experiments considered - DES [37], LSST [38], the SuperNova/ Acceleration 
Probe (SNAP) [39], and Planck [40]. We present our results in Sec. IV along with a discussion of their dependence 
on the assumptions made in the analysis. We conclude with a summary in Sec. V. 
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II. PARAMETRIZING THE MODIFICATIONS OF GRAVITY 



We study the evolution of linear matter and metric perturbations in a general metric theory of gravity. Assuming 
that the background evolution is correctly described by the flat Friedmann- Roberts on- Walker (FRW) metric, we focus 
on scalar perturbations and work in the conformal Newtonian gauge, so that the perturbed line element is given by 

ds^ = -a\rj)[{l + 2^'{x,j]))dr]^ - {l-2<i>{x,ri))dx^] , (1) 

where rj is the conformal time. In the rest of the paper, all the perturbed quantities are presented in Fourier space, 
dots indicate derivatives with respect to r], and Ti. = d/a (as opposed to H, which contains a derivative with respect 
to cosmic time t, H = a^^da/dt). We use the standard notation for the energy momentum tensor of the matter fields, 
which to first order in the perturbations, assumes the following form 

T% = -p{l+S), 

T", = -(/> + (2) 
= (P + 5P)<5}+7r;, (3) 

where 5 = Sp/p is the density contrast, v the velocity field, SP the pressure perturbation and tt* denotes the tracelcss 

component of the energy-momentum tensor. Finally, we define the anisotropic stress a via {p+P)a = —{k kj — ^<5p7r* . 

In GR, the linearized Einstein equations provide two independent equations relating the metric potentials and 
matter perturbations, the Poisson and anisotropy equations, respectively: 

^'* = -:^pA, (4) 



2Mf, 



2 



k'{'P-^) = ^{p + P)<j, (5) 

where pA = p5+3^{p+P) V is the comoving density perturbation. In the ACDM and minimally coupled quintessence 
models, the anisotropic stress is negligible at times relevant for structure formation, and we have = $. 

In models of modified gravity, as well as in more exotic models of dark energy, the relation between the two New- 
tonian potentials, and between the potentials and matter perturbations, can be different [19, 41, 42]. We parametrize 
the changes to the Poisson and the anisotropy equations as follows: 

^'* = -^/^(«,fc)pA (6) 

|-7(a,fc), (7) 

where p{a, k) and 7(0, k) are two time- and scale-dependent functions encoding the modifications of gravity and/or 
the contribution of an exotic dark energy fluid. Note that we have chosen to define p via the Poisson equation (6) 
written in terms of 4", the perturbation to the time-time component of the metric. This choice is natural, as it is 
that enters the evolution equation for CDM density perturbations on sub- horizon scales^: 

S + nS + k^"^ ^0 . (8) 

The gravitational potentials are not going to be observed directly. In order to see the effect of the modifications 
on observable quantities, such as Galaxy Counts (GC), Weak Lensing shear (WL) and CMB, we modify the publicly 
available Code for Anisotropies in the Microwave Background (CAMB) [43, 44]. This requires implementing Eqs. (6) 
and (7) in the synchronous gauge used in CAMB. The details of this procedure are given in Appendix A, where we 
also demonstrate that our method respects the super-horizon consistency condition [42, 45]. 

In addition to the effects of modified gravity, the evolution of all cosmological perturbations depends on the back- 
ground expansion. We restrict ourselves to background histories consistent with the fiat ACDM model. One reason 
for this choice is that ACDM is currently the best fit to available data. Another reason comes from the fact that in 
the recently popular models of modified gravity, e.g. the scalar-tensor models, f{R), nDGP and Degravitation (but 



^ Our rescaling of the Newton's constant fi corresponds to G^/G in the notation of [32] - their <I> is our 'if, and vice versa. 
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not sDGP), the expansion history is efFeetively the same as in ACDM, and the main differences between models arise 
from the evolution of cosmic structure. Finally, by restricting ourselves to ACDM-like expansion histories, we can 
better distinguish the effects of modified growth. Wc do not, however, fix the values of the cosmological parameters 
flmh^ , ^bh^ and /i, nor the spectral index Ug or the optical depth r. These cosmological parameters will be varied 
along with the modified growth parameters. In principle, it would be interesting to also allow for a non-zero spatial 
curvature and for variations in the effective dark energy equation of state. We have included these effects in our 
forecasts for ^{a, k) and 7(0, k) based on the PCA approach [35]. 

To arrive at a suitable parametrization of the functions /^(a, k) and 7(0, fc), we note that models of modified gravity 
typically introduce a transition scale which separates regimes where gravity behaves differently. For example, in f{R) 
and scalar-tensor models, the functions /x and 7 are equal to unity at early times and on large scales {i.e. on scales 
that are larger than the characteristic scale of the model) , and they transition to a modified value on smaller scales 
and late times. We can mimic this time- and scale-dependent transition via the following functions of the variable 
k'^a" 



lJ.{a,k) 
7(a, k) 



1 + k^a' 
1 + A? k'^a" 
1 + 132X1 k^a' 



(9) 
(10) 



where the parameters Xj have dimensions of length squared, while the pi represent dimensionless couplings. The 
expressions (9) and (10) coincide with the scale-dependent parametrization introduced in [32]. It is easy to show 
that this parametrization follows from scalar-tensor theories, of which the f{R) models are an example. In these 
models, one makes a distinction between the so-called Jordan frame, with the metric 5^1/, (where matter falls along 
the geodesies and the action for gravity is modified), and the Einstein frame, with the metric g^i^, in which the Einstein 
form of the action is preserved but there is an additional scalar field non- minimally coupled to gravity [46-48]. In 
what follows, we adopt the usual convention of indicating Einstein frame quantities with a tilde. The two frames are 
related through a conformal mapping: 



(11) 



where ai(<j)) is defined in Eq. (12), and the index i indicates different matter components. In the Einstein frame, the 
action is a Chameleon- type [9] action for the massive field cj) coupled to matter fields via couplings 0,(0) 



_£i?__gM.(V^,/,)V,0-y(0) 



(12) 



For f{R) theories, the coupling is universal and linear, a ~ ■\/2/3 however, in a general scalar-tensor theory, the 
coupling(s) can be a non-linear function(s) of the field (j). In what follows, we will focus on the growth of structure 
and consider primarily cold dark matter, working with a single coupling a{(j)). 

The behavior of cosmological perturbations in general models of coupled dark energy has been studied in detail 
in [49-52]. The scalar field does not contribute any anisotropic stress, neither do matter fields, so the Newtonian 
potentials in the Einstein frame are equal: $ = However, the effective potential acting on dark matter particles 
has an extra contribution due to the interaction of matter with the field [50]. Neglecting baryons and radiation, 
one obtains the following Poisson equation 



eff 



3 0?Pn 
"2~fc2" 



1 



\a'^Y{k) 



(13) 



where Y{k) = k"^ /{k"^ + a^m?) is the Yukawa term, m is the time-dependent effective mass of the scalar field, and 
primes denote differentiation w.r.t. the field 0. Following the mapping prescription described in Appendix B, we can 
now map these equations back to the Jordan frame. As a result of this mapping, we obtain the following effective 
Newton constant p and ratio between the potentials <& and 5", 7: 



1 



/i(a, k) 



lia,k) = 



1 



2" ) 



2 I a^m 



— 2" 



2 / a^m^ 



fe2 



(14) 
(15) 
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In writing the expression for fj,, we neglected an overall pre- factor of e~'^°'^'^\ which corresponds to a time-dependent 
rescaling of the Newton constant. In the specific case of f{R) theories it corresponds to 1 -I- fa, where = df /dR. 
The consistency of scalar-tensor theories with locaP and cosmological tests requires that e"'^"^'^' ~ 1, with departures 
from unity being outside the reach of current and upcoming cosmological probes. Therefore, it is safe to approximate 
it with unity. 

The expressions (14) - (15) are equivalent to the parametrization (9) - (10) once the parameters {s,Xf,/3i} are 
chosen as follows 

^(1+^2) ^ ^ 

m 





1 

ml 




xl = 


1 

ml 


H 


/3i = 


1 + 


a 

'Y 


/32 = 


2- 




2 + 


a'2 ■ 



(16) 

In deriving Eqs. (16), we have assumed a' ~ const. This is exact in /(i?) theories, where the coupling a(0) is a linear 
function of the field. In more general scalar-tensor theories, the coupling can be a non-linear function of the scalar 
field. However, the value of the field (/), typically, does not change significantly on the time-scales associated with the 
epoch of structure formation. It is, therefore, reasonable to approximate a' with a constant for our purposes. 
From Eqs. (16), it is clear that in the case of scalar-tensor theories, the parameters {A?,/3i} are related by 

/3i = | = 2-/?,|. (17) 

By fitting forms (9) and (10) to data and checking the validity of the relations (17), one can test and potentially rule 
out a large class of scalar-tensor models of cosmic acceleration. 

In our forecasts, we will need to assume particular fiducial values of {s,Af,/?i}. One type of fiducial cases we 
consider are based on f{R) theories. The other type is motivated by the Chameleon scenario. In any scalar-tensor 
theory, the parameter s is determined by the time-evolution of the mass of the scalar. The non-minimal coupling 
leads to the dependence of the effective mass of the scalar on the local energy density of non-relativistic matter. The 
effective potential for the scalar field is 



which gives an effective mass 



m" 



Feff(</)) =F(0)H-p„,e«"('^) , (18) 



V;'4('/'min) = V"-n (a" + a'') V , (19) 



calculated at the minimum 0min of the potential V^ff. Following [56], the time dependence of the mass (19) can be 
approximated as 



m _ 1 V" ■ 

^ o T777 Vrain ■ 

m 2 V 



(20) 



For tracking- type potentials described by an inverse power-law, i.e. V ^ (j) "j the mass evolves as 

m ~ a-3(«+2)/2("+i) (21) 

corresponding^ to 1 < s < 4. Cosmologically viable f{R) models correspond to s ~ 4. The scalar degree of freedom 
there is represented by the function and its mass is a time dependent function set by « f^]^ [14, 16]. From 



^ Even if f(R) models can be made to satisfy solar system constraints, they may still lead to singularities when one considers formation 

of compact objects, such as neutron stars [53-55] . Avoiding these singularities would require additional fine-tuning. 
* It is easy to show that s must be in this range for any V{4>) ~ 0", for both positive and negative a. 
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numerical simulations [16], we find that fjm oc a^, corresponding to s ~ 4. This can be understood analytically using 
the following argument. A viable /(i?) must be a fairly slowly- varying function of R, and in the relevant range of R 
we can always approximate it with a power law: /(i?) ^^01+ a2i?" up to the leading order in R. This function gives 
fjifj ~ n{n — l)f /E? ~ const • R~^. The mass of the field is calculated at the minimum of the scalaron potential, 
where R = K?p, therefore, we have /ij^ cx oc a^, which is what we have found numerically. In the case of 
tracking-type quintessence potentials, usually considered in Chameleon scenarios, one needs small integer values of 
the index n > 0{1) in order to have the appropriate mass scale [56]. Therefore, s ^ 2 can be seen as a typical value 
for Chameleon models. 

The fiducial values of Pi and /32 can be derived from a' using Eqs. (16). f{R) models correspond to a linear coupling 
with a' = ■\/2/3. In this case, Eqs. (16) give 

/3i = I /32 = i . (22) 

In Chameleon models, a' is a free parameter that is typically assumed to be ~ 0(1). 

Finally, the mass scale today, mo, is a free parameter in all scalar-tensor models. In the f{R) case, it is constrained 
from below by requirements of consistency with local and cosmological tests. In particular, local tests of gravity [14] 
set a lower bound, mo > 10~^ Mpc~^, which corresponds to < 10^ Mpc^. This bound can be relaxed to mo > 10"'^ 
Mpc~^, corresponding to < 10^ Mpc^, if we consider only cosmological tests, neglecting local constraints [57]. In 
the case of Chameleons, the constraints are somewhat weaker due to the additional freedom in choosing a' . However, 
if a' ^ 0(1), then the bounds on toq should be comparable to those in the f{R) case. 

In our forecasts, we use two f{R) fiducial cases, corresponding to two choices of the mass scale toq- Namely, in 
both f{R) cases, we take s = 4, Pi = 4/3 and P2 = 1/2, and try two different sets of values for (A^, X^)- Model I with 
A2 = PiXi = 10^ Mpc^ and Model II with A| = /3iA^ — 10^ Mpc^. Wc also consider two fiducial models corresponding 
to the Chameleon case, both with s ~ 2, Pi = 9/8 and P2 = 7/9, which correspond to a' = 0.5. We then use the 
same fiducial values for [XfjX^) as in the f{R) case. We refer to these two models as Model III and Model IV. In 
practice, rather that varying dimensionful parameters A|, we work with log(A|/Mpc^). To summarize, we consider 
the following fiducial cases: 

• Model I, with s = 4, /3i = 4/3, P2 = 1/2, log(Ai/Mpc^) = log(/3iAf /Mpc^) = 3; 

• Model II, with s = 4, /3i = 4/3, P2 = 1/2, log(A2/Mpc2) = log(/3iAf /Mpc^) ^ 4; 

• Model III, with s = 2, Pi= 9/8, p2 = 7/9, log(A^/Mpc^) = log(/3iA?/Mpc^) = 3; 

• Model IV, with s = 2, Pi= 9/8, P2 = 7/9, log(A2/Mpc2) = log(/3iA^/Mpc2) = 4. 

In Fig. 1, we plot the functions fi{a, k) and 7(0, k) corresponding to the four fiducial models considered in this work. 
In the ACDM case, both functions are constant and equal to unity. For Models I-IV, the two functions approach 
unity at small values of fc^a**, while at larger k'^a", they transition to values determined by Pi and P2- 

In the recent literature, there has been a growing interest in parametrizing modifications of gravity. The parametriza- 
tion we have used is equivalent to the scale-dependent form proposed in [32]. In the same paper, the authors also 
consider a scale-independent parametrization, in which fi and 7 arc not independent functions. A scale-independent 
parametrization of the ratio of the potentials was also introduced in [58], although with a fixed time-dependence. 
There is, however, no physical model of modified gravity which corresponds to a scale-independent modification of 
growth. Indeed, to avoid serious confiicts with Big Bang Nucleosynthesis (BBN) and CMB data, the modified gravity 
theories must reduce to GR at early times. They should, however, introduce novel physics at late times on cosmologi- 
cal scales, while obeying local constraints of gravity on non-linear scales. These requirements introduce characteristic 
transition scales into the theory that lead to a peculiar scale-dependent behavior of perturbations. A parametriza- 
tion describing all these three regimes, including a possible modified time-dependence on super-horizon scales, is the 
Parametrized Post-Friedmann Framework (PPF) of [29]. It consists of three functions of time and space, plus one 
parameter, plus a parametrization of the non-linear regime. In this work, we are interested in modifications of gravity 
that primarily become important on sub-horizon cosmological scales. For this purpose, our parametrization, which 
has five parameters, is sufficient and more economical. That is, our parametrization, by design, reduces to GR on 
super-horizon scales and evolves into a modified version of gravity only below a certain sub-horizon scale, such as the 
Compton wavelength of the scalar field. Furthermore, we do not consider the non-linear scales, because describing 
them properly requires input from A^-body simulations. This can be done in specific models, such as particular 
types of f{R) [59], but would be ill-defined if one tried to model them based on a general parametrization of linear 
gravitational potentials, such as Eqs. (9) and (10). 
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FIG. 1: The rescaling of the Newton constant ^{a,k) and the ratio of Newtonian potentials 7(0, fc), plotted as a function of 
the redshift for the four fiducial models used in our Fisher analysis. The contours are lines of constant ^ and 7. The first two 
models correspond to f{R) fiducial cases, with s = 4, the coupling /?i = 4/3 and the mass scale of \l = lO^Mpc^ for model I 
and A| = lO^Mpc^ for model IL Models III and IV correspond to Chameleon theories with s = 2, the coupling /3i =9/8 and 
the mass scale of A2 — lO'^Mpc^ and A2 — lO'^Mpc^, respectively. 



The parametrization we use gives an accurate description of the effects of scalar-tensor theories (and in gen- 
eral models of coupled quintessence) on linear scales. Describing models such as DGP would require the introduc- 
tion of time-dependent modifications on super-horizon scales (and therefore the PPF framework would be a better 
parametrization in that case). However, we expect that Eqs. (9) and (10) can capture some of the features of DGP 
if applied to transitions on larger scales. In [30], the authors have proposed a framework to test GR on cosmological 
scales; the parametrization they use is suitable for braneworld theories, while less suited for describing models of the 
/(i?) type. As already mentioned in the introduction, a very promising model-independent procedure to test the 
relation between the metric potentials which uses the cross-correlation between peculiar velocities and galaxy counts, 
and galaxy counts with Icnsing, was introduced in [28]. That method probes 7 but not /x, and extracting peculiar 
velocities is notoriously difficult. Finally, in a recent paper [52], the authors considered constraints from CMB and 
LSS on general theories of interacting dark matter and dark energy. 



III. TOMOGRAPHIC OBSERVABLES 



In the ACDM model, the sub-horizon evolution of gravitational potentials and the matter density fluctuations are 
described by a single function of time- the scale-independent growth factor g{a). Namely, for all Fourier modes 
that have entered the horizon by some epoch a^, we have A(a, A;)/A(ai, /c) = 17(0) and the potentials evolve as 
^{a,k)/'^{ai,k) = ^{a, k)/^{ai, k) = g{a)/a. In models of modified gravity, on the other hand, the dynamics of 
perturbations can be richer and, generically, the evolution of \t' and A can be described by different functions of 
scale and time. By combining different types of measurements, one can try to reconstruct these functions, or at least 
put a limit on how different they can be. In what follows, we give a brief introduction into the relation between the 
different types of observables and the gravitational potentials they probe. For a more thorough review of the various 
ways of looking for modifications in the growth of perturbations we refer the reader to [60] . 

Galaxy Counts (GC) probe the distribution and growth of matter inhomogcneitics. However, to extract the matter 
power spectrum, one needs to account for the bias, which typically depends on the type of galaxies and can be both 
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time- and scale-dependent. On large scales, where non-linear effects are unimportant, one can use a scale-independent 
bias factor to relate galaxy counts to the total matter distribution. This relation becomes increasingly complicated and 
scale-dependent as one considers smaller and smaller scales. In principle, the bias parameters can be determined from 
higher order correlation functions [61-63]. In this work we stick to linear scales where the bias is scale-independent. On 
sub- horizon linear scales, the evolution of the matter density contrast is determined by Eq. (8). Hence, measurements 
of GC over multiple redshifts can provide an estimate of 'J as a function of space and time, up to a bias factor. A more 
direct probe of the potential would be a measurement of peculiar velocities, which follow the gradients of 5*. Such 
measurements would be independent of uncertainties associated with modeling the bias. Peculiar velocity surveys 
typically use redshift-independent distance indicators to separate the Hubble flow from the local flow, and nearby 
SNela are therefore good candidates; a number of surveys, like the 6dFGS [64] and the 2MRS [65], use galaxies. An 
interesting alternative is offered by the kinetic Sunyaev-Zel'dovich effect of clusters [66], that arise from the inverse 
Compton scattering of CMB photons off high-energy electrons in the clusters. This effect provides a useful way of 
measuring the bulk motion of electrons in clusters, hence the peculiar velocity of clusters, but it is limited by low 
signal-to-noise ratio. Current measurements of peculiar velocities are limited in accuracy, and at this point it is not 
clear how to forecast the accuracy of future observations. Therefore we did not include them in our observables, even 
though they are a potentially powerful probe [67]. 

In contrast to galaxy counts and peculiar velocities, which respond to one of the metric potentials, namely \E', 
Weak Lensing (WL) of distant light sources by intervening structure is determined by spatial gradients of ($ -I- ^'). 
Hence, measurements of the weak lensing shear distribution over multiple redshift bins can provide an estimate of 
the space and time variation of the sum of the two potentials. In the ACDM and minimally coupled models of dark 
energy, the two metric potentials coincide and therefore WL probes essentially the same growth function that controls 
the evolution of galaxy clustering and peculiar velocities. In models of modified gravity, however, there could be a 
difference between the potentials, corresponding to an effective shear component, also called a "gravitational slip". 

Measurements of the Integrated Sachs- Wolfe effect (ISW) in the CMB probe the time dependence of the sum of 
the potentials: <& + '!'. The ISW effect contributes a gain (loss) of energy to CMB photons traveling through decaying 
(growing) potential wells. It contributes to the CMB anisotropy on the largest scales and is usually extracted by 
cross-correlation of the CMB with galaxies [68-70] . 

By combining multiple redshift information on GC, WL and CMB, and their cross-correlations, one can constrain 
the differences between the metric potentials and the space-time variation of the effective Newton constant defined 
in the previous section. Ideally, the experiments would provide all possible cross-correlations, between all possible 
pairs of observables, in order to maximize the amount of information available to us. In practice, however, it can be 
difficult to obtain these cross-correlations, since their measurements require that each of the individual fields (CMB, 
GC, WL) be measured on the same patch of sky. This will be addressed with near and distant future tomographic 
large scale structure surveys (like DES [37], PAN-STARR [71] and LSST [38]). We have made separate forecasts using 
various combinations of data sets, and all possible cross-correlations are considered (CMB-WL, CMB-LSS, LSS-WL). 
In this manner, we can determine how well modified gravity will be constrained with varying degrees of experimental 
difficulty and sophistication. 



Consider a two-point correlation function C^^ [9) = C^^ {\ni — h2\) = {X{hi)Y{h2)) between two 2D-fields, X 
and Y, measured by an observer looking at the sky. Here n is a direction on the sky and cos 6' = fii • n2. C^^ {&) can 
be expanded in a Legendre series: 



where the first two terms in the series (the monopole and dipole contributions) are coordinate dependent and should 
vanish in the CMB frame of a homogeneous and isotropic universe. The expansion coefficients Cf^ can be expressed 
in terms of the primordial curvature power spectrum and the angular transfer functions lf'^(k) as 



A. Angular spectra 




(23) 




(24) 



where 




(25) 
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FIG. 2: A schematic representation of all the 2-point correlations used in this work. M is the number of redshift bins used 
for galaxy counts (GC), A'^ is the number of bins used for weak lensing (WL), T and E stand for temperature and E-mode 
polarization of the CMB. The white blocks corresponding to the cross-correlations between E and Gi, and between E and WL; 
indicate that those cross-correlations were not considered in this work. 



and similarly for . A detailed derivation of the above expressions is given in Appendix C. The observable quantities 
for which we need to evaluate the I('s are the galaxy distributions at several ranges in redshift (the tomographic redshift 
bins), the maps of lensing shear for different bins, and the CMB temperature anisotropy. 

The distribution of galaxies is expected to trace the distribution of dark matter up to a bias factor which quantifies 
selection effects specific to the type and color of the galaxies. The bias factor becomes increasingly scale-dependent on 
smaller scales and its modeling, e.g. using the halo model [72, 73], involves assigning on the order of five parameters 
to each photometric bin. The halo model assumes the validity of General Relativity and would not be directly 
applicable for testing modified gravity. Also, the modifications of growth that we are trying to detect would be 
highly degenerate with the variations of the bias parameters. Higher order correlation functions can, in principle, 
significantly help in reducing this degeneracy [61-63]. Still, this would require calculating three-point statistics of 
cosmological perturbations in modified growth models, which is beyond the scope of this paper. Instead, we work 
under the assumption that on large scales, the bias can be treated as scale-independent and can be modeled with one 
free parameter bi for each redshift bin i. With this assumption, the corresponding angular transfer functions can be 
expressed in terms of the dark matter density contrast as 

/f (fc) = hcsn I * dzWGAz)ji[kr{z)\5{k,z) , (26) 
Jo 

where Wctiz) is the normalized selection function for the «th redshift bin, S{k,z) is the density contrast transfer 
function, and csn = —9/10 for the adiabatic initial conditions [74, 75] that are assumed throughout this work. 
For weak lensing, the relevant le's are given by 

= c*TC / * dzW^Mn[kr{z)]{^ + c$*<l) , (27) 
Jo 

where c-qj-n = 3/5, c<E,* = (1 + 2i?^/5), where = (7iV^/8)(4/ll)''/3 [74], is the number of flavors of relativistic 
neutrinos, and W^^ (z) is the window function for the ith bin of sheared galaxies with a normalized redshift distribution 

W,^ (z) = r (^') . (28) 

Jz nz) 
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FIG. 3: Some representative growth functions (time evolution divided by the corresponding initial values, set at z = 30) that 
contribute to the integration kernels of the Ce's for ISW, WL, GC, and their cross-correlations. The corresponding Ct's are 
shown in Fig. 4. The solid line corresponds to GR and the other lines are for Model II (an f{R) fiducial model described in 
Sec. II). The different types of line types/colors represent four different k modes and are explained in the legend. 



The transfer functions for the CMB temperature anisotropy receive contributions from the last-scattering surface (at 
z ~ 1100) and from more recent redshifts via the Integrated Sachs- Wolfe (ISW) effect. By design, our modifications 
of gravity should be negligible at recombination. Therefore, their only imprint on the CMB will be via the ISW effect. 
For the ISW contribution to the CMB. wc have 



rISW 



(k) = 



dzt 



(29) 



where r(z) is the opaqueness function. 

We numerically evaluate the transfer functions to obtain using MGCAMB: a code developed by us'' for 

studying modified growth based on CAMB [43]. The details of the implementation of the parametrization (6) in 
CAMB, which involved conversion to the synchronous gauge, are given in Appendix A. To save computing time, we 
used the fact that modifications of GR are negligible at z > 30, and until then the usual CAMB solver was used. 
From z = 30 and on we continue by solving Eqs. (A13— A14). 

As illustrated in Fig. 2, a joint analysis of CMB and data from a tomographic lensing survey with AI GC redshift 
bins and TV WL bins can give us a total of 3 + M{M +l)/2 + N{N +l)/2 + M + N + MN different types of C/s from 
CMB, GC, WL, GCxCMB, WLxCMB and WLxGC, respectively (we do not correlate CMB polarization with GC 
and WL). For example, combining Planck with DES, with M = 4 GC bins and N = 4 WL bins, gives us 47 different 
types of spectra. A combination of Planck with LSST, with 10 GC bins and 6 WL bins, gives us 155 different Cg's. 



5 Publicly available at http://www.sfu.ca/~gza5/MGCAMB.htmI 
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FIG. 4: Dimensionless power spectra l[t+l)C^^ /2it, where X,Y = {T, GC, WL} for GR (solid line) and for our Models I-IV 
(line types explained in the legend) , for Planck and a few representative choices of GC and WL redshift bin pairs expected from 
LSST. On the plots, stands for the i-th WL bin, and gi stands for the z-th GC bin. The assumed redshift distribution of WL 
and GC sources for LSST, along with their partition into photometric bins, is shown in Fig. 5. For each group of spectra, we 
pick one with the largest deviation from the GR prediction and plot its relative difference w.r.t GR in the lower part of each 
panel. The straight dashed lines indicate the approximate scale at which non-linear corrections become significant. We have 
only used the parts of the spectra that can be accurately described by linear theory. 



For a given function / ~ f{k, z) where k and z arc conioving wavenumbcr and redshift respectively, we can define 
a corresponding growth function D as 

In Fig. 3, we plot these growth functions corresponding to the kernels of several representative C^'s for ISW, WL, GC 
and their cross-correlations. In GR, all of these quantities grow in a scale-independent manner, while in the modified 
scenarios considered here their growth is enhanced in a scale-dependent way. We see that the enhancement is more 
pronounced on smaller scales (larger fc), as expected in scalar-tensor scenarios [16]. 

In Fig. 4, we show some of the spectra expected from a combination of Planck and LSST for the four fiducial models 
considered in this paper. Shown are the auto- and cross-correlated spectra for two representative WL and GC redshift 
bins, with central values z{ei) « 0.3, ^(£2) « 0.8, 2(53) « 0.4, 2(35) « 0.9, where et stands for the i-th WL bin, and gi 
stands for the i-th GC bin. The full redshift distribution of WL and GC sources we assumed for LSST, along with 
their partition into photometric bins, is shown in Fig. 5. 

We show only the parts of the spectra that correspond to the linear cosmological regime. Including higher or 
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smaller scales, would require us to account for non-linear effects which, strictly speaking, is not allowed within our 
framework. To accurately model growth on non-linear scales, one needs input from N-body simulations, which can only 
be performed for specific modified gravity theories. The fact that we are not testing a specific model, but constraining 
a general departure from GR, defined in terms of linear perturbation variables, precludes us from having a reliable 
description of non- linear corrections. Simply applying the analytical corrections developed under the assumption of 
GR, such as by Peacock and Dodds [76] and Smith et al [77], can lead to significant errors [59]. For this reason, 
we restrict ourselves to the linear regime by cutting off the spectra at /max ~ 0.2 hx{zs), where xi^s) is the 

comoving distance to the redshift bin at z = Zj. This cutoff roughly corresponds to fc ~ 0.2 /iMpc~^ at z = 0. There 
is certainly a wealth of information about MG parameters on smaller scales. For comparison, using fc ~ 0.1 /iMpc~^ 
at z = as a cutoff degrades the LSST constraints on MG parameters by a factor of ~ 2. However, while it would be 
tempting to include information from even smaller scales, as mentioned above, it would make predictions obtained 
using linear theory unreliable. 



B. Fisher matrices 



In order to determine how well the surveys will be able to constrain our model parameters, we employ the standard 
Fisher matrix technique [78] . The inverse of the Fisher matrix Fab provides a lower bound on the covariancc matrix of 
the model parameters via the Cramer- Rao inequality, C > F^^ [78]. For zero-mean Gaussian-distributed observables, 
like the angular correlations Cf^ introduced in Section III A. the Fisher matrix is given by 



= E^-f^C-^C-), (31, 



where Pa is the a^^ parameter of our model and Cg is the "observed" covariance matrix with elements Cf^ that 
include contributions from noise: 

ATf^ . (32) 

The expression (31) assumes that all fields X(n) are measured over contiguous regions covering a fraction /sky of the 
sky. The value of the lowest multipolc can be approximately inferred from £,„in ~ 7r/(2/sky). It is also possible to 
write expressions for separate contributions to the Fisher matrix from particular subsets of observables. For example, 
for angular spectra ' ^, corresponding to N fixed pairs of fields {Xi,Xj) we can write 

^If = E (2^+1) E E ^^[^^"]'^^^%^ ' (33) 

fcf„,„ {.y} = l {,„„} = ! 

where the covariance matrix Cl"*^ has elements [72] 

j^sub]{»i}{mn} ^ ^X,X^^X,X,, ^ ^X^X^^X.X^ ^g^^ 

Eqs. (31) and (33) become the same in the limit of summing over all possible pairs {Xi, Xj). 

In principle, the noise matrix includes the expected systematic errors. Systematics are, however, notoriously 

difficult to predict, and are often ignored in parameter constraint forecasts. For both surveys considered in this work, 
DES and LSST, achieving their science goals requires bringing the systematic errors under the noise level [37, 79]. 
Therefore, we take Nf-^ to represent statistical noise only. Assuming uncorrelated Poisson noise on the galaxy 
overdensity in each bin (Gi) and shear fields (k^), the noise is given by [72] 

'h 

Nf'"' = 0, (35) 

where ^rms is the expected root mean square shear of the galaxies, and Uj is the number of galaxies per steradian in 
the jiYi redshift bin. 
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For CMB. the temperature anisotropy T and the E-mode polarization are measured in several frequency channels. 
When combined, the total noise is less than that in any of the channels 



N. 



YY 



-1 -1 

-1 



(36) 



where is the noise on a measurement in a given channel and Y denotes either T or E. Assuming the CMB 

experiment sees Gaussian beams, this autocorrelation noise is given by 

j^YY ^ f ^I^FWHUA gf(«+l)e|„HM,c/81n2 ^ ^37^ 
V TcMB J 

where is the standard deviation in observable Y (either T or E) in channel c, and 6'fwhm,c is the FWHM of the 
beam, in arc-minutes, at the given frequency channel [80] . 

For supernovae, the direct observable is their redshift-dependent magnitude 

m(z) = 7W + 5 log + 25 (38) 

where Ai is the intrinsic supernova magnitude and the luminosity distance (in Mpc) defined as 

rfL(.) = (l + ^)^^;^ , (39) 

where H{z) is the Hubble parameter with a current value of Hq. The information matrix for SNe observations is 

pSN _ 1 dm{z^) dm{zi) 

Z^a{z,r dpa dp, ■ ^ ' 

where the summation is over the redshift bins and am{zi) is the value given by Eq. (45) at the midpoint of the i-th 
bin. 

Thus, given a set of theoretical covariance matrices over a given multipole range, and the specifications for the ex- 
pected noise in particular experiments, we can compute the Fisher matrix. The derivatives with respect to parameters 
are computed using finite differences. The values of the finite differences are given in Appendix D. Upon inverting the 
Fisher matrix, we find the theoretical lower limit on the covariance of the parameters of the model. In this manner, 
we can forecast how tightly upcoming experiments will be able to constrain modified growth. 



C. Experiments 



The data considered in our forecasts include CMB temperature and polarization (T and E), weak lensing of dis- 
tant galaxies (WL), galaxy number counts (GC), their cross-correlations, and SNe observations. We assume CMB T 
and E data from the Planck satellite [40], the galaxy catalogues and WL data by the Dark Energy Survey [37] 
and Large Synoptic Survey Telescope (LSST) [38], complemented by a futuristic SNe data set provided by the 
ongoing Canada-France-Hawaii Telescope's (CFHT) Supernovae Legacy Survey (SNLS) [81, 82], the Nearby Su- 
pernovae Factory (NSNF) [83, 84], and a future Joint Dark Energy Mission (JDEM) space mission, such as the 
Supernovae/ Acceleration Probe (SNAP) [39]. 



1. DBS and LSST 



To describe the expected redshift distribution of galaxies, we take the total galaxy number density to be given by 

Naiz) cx z\xp{-z/zof , (41) 

which is a a slight modification of the model due to Wittman et al. [85]. This function represents a compromise 
between two trends; more distant galaxies are harder to see, so they are less likely to be a part of the field we observe, 
but as one goes to higher redshifts, there is more volume of space available and thus one expects there to be more 
galaxies. The parameter zg depends on the experiment and defines the redshift at which the most galaxies will be 
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FIG. 5: The assumed redshift distribution and the partition into photometric bins for the galaxy counts (GC) and weak lensing 
(WL) for DES and LSST. The bias parameter for each GC bin is shown with a thin red solid line in the lower panel. 



observed. The value of Zq for a given experiment depends on how faint an object it can observe. The window function 
is suitably normalized so that the total number of galaxies matches that expected from a given survey. The galaxies 
can be divided into photometric redshift bins, labelled with index i, 

7Vg(z) = ^iVG.(z). (42) 

i 

In our analysis, we assume that the photometric redshift errors are Gaussian distributed, and that their rms fluctua- 
tions increase with redshift as cr(z) = CTmax(l + + Zmax)- The bin sizes are chosen to increase proportionally to 
the errors. The resulting photometric redshift distributions are given by. 



f Zi-l - z\ f Zi~ z 

erf c — ;= — ertc —= 

\V2<j{z)J \V2a{z) 



(43) 



where erfc is the complementary error function. For a given photometric redshift bin, the normalized selection function 
that appears in Eq. (26) is given by 

WaAz) = ^ (44) 

where A^* is the total number of galaxies in the z-th bin. 

DES is a project aimed at studying the nature of the cosmic acceleration, and is planned to start observations 
in September, 2009 [37]. DES includes a 5000 square degree multi-band, optical survey probing the redshift range 
0.1 < z < 1.3 with a median redshift of z = 0.7 and an approximate l-cr error of 0.05 in photometric redshift. In our 
simulation, for both WL and galaxy counts, we assume a sky fraction /sky = 0.13, and an angular density of galaxies 
Ng = 10 gal/arcmin^. We also assume 7rms = 0.18 -I- 0.042 z, which is the rms shear stemming from the intrinsic 
eflipticity of the galaxies and measurement noise, and the photometric redshift uncertainty given by a{z) = 0.05(l + z). 

LSST is a proposed large aperture, ground-based, wide field survey telescope [38]. It is expected to cover up to half 
of the sky and catalogue several billion galaxies out to redshift z ~ 3. For LSST forecasts, we adopt parameters from 
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TABLE I: The relevant parameters for Planck [88]. We use the three lowest frequency channels from the HFI. 
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TABLE IL The binned redshift distribution of type la supernovae used in our analysis. The redshifts given are the upper limits 
of each bin. Magnitude errors cFm[z) are evaluated at bin midpoints. 



the recent review paper by the LSST coUaboration [79]. Namely, we use /sky = 0.5, Nq ~ 50 gal/arcmin^ for both 
WL and counts, = 0.18 + 0.042 z, and a{z) = 0.03 (1 + z). 

For both DES and LSST, we take the GC photometric bins to be separated by ha{z). This leads to four redshift 
bins for DES and ten for LSST. For WL (source) galaxies, we use four bins for DES and six for LSST. In principle, we 
could have used a larger number of WL bins for LSST, and it would be interesting to investigate how the constraints on 
modified growth improve with finer binning. It is known, for example, that finer binning does not improve constraints 
on the dark energy equation of state w{z) [72]. The assumed galaxy number distributions for both experiments are 
shown in Fig. 5. 

The photometric redshift errors used in our analysis should be seen as the optimistic values for the respective 
experiments. For a discussion of potential sources of systematic errors in photo-z estimates, the reader is referred 
to [86, 87]. We assume that the systematic errors for both DES and LSST will be comparable or smaller than the 
statistical errors. While achieving this level of systematics control will be a huge challenge, the white papers for both 
experiments [37, 79] describe it as possible and, in fact, necessary for accomplishing the science goals set by these 
experiments. 

2. Planck 

The Planck mission [40] of the European Space Agency (ESA) is currently expected to launch in the Spring of 2009. 
Planck will image the full sky with a sensitivity of AT/T ~ 2 x 10~^, angular resolution to 5', and frequency coverage 
of 30-857 GHz [88]. The angular resolution of Planck will be three times superior to that of NASA's WMAP satellite, 
and the noise lowered by an order of magnitude at around 100 GHz. These significant improvements will permit more 
accurate measurements of the CMB temperature and polarization power spectra, allowing for a better determination 
of the cosmological parameters. We use the expected sensitivity and resolution parameters for the lowest three Planck 
HFI channels based on a 14 month mission [88]. We list the parameters relevant to our Fisher analysis in Table I. 

3. Supernovae and the parameter priors 

In our forecasts, we assume spatially flat geometry and expansion histories, consistent with the effective equation 
state of dark energy equal to —1. In addition to the five modified growth parameters, and the M bias parameters, 
we vary the Hubble constant h, cold dark matter density flch^, the baryon density Vt^jh?, the optical depth t, the 
scalar spectral index Ug, and the amplitude of scalar perturbations Ag. Their fiducial values are taken to be those 
from the WMAP 5-year data best fit [89]: n^h? = 0.023, f^^/i^ = 0.11, = 0.72, r = 0.087, = 0.963. The fiducial 
values for bias parameters are motivated by the parametrized halo model described in [72] , and we show them in the 
lower panel of Fig. 5. Imposing a prior on the value of h from the Hubble Space Telescope (HST) did not make a 
noticeable difference in our results. 

In addition, to better constrain the background expansion parameters, we include data from a futuristic set of SNe 
luminosity distances from SNLS, NSNF and SNAP. SNe observations determine the magnitude of the SNe, m(z). 
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FIG. 6: The 68% confidence contours for our five modified growth (MG) parameters for Model I, based on f(R) (described 
in Sec. II), as constrained by several difi'erent combinations of correlation functions from various experiments (shown in the 
legend). As described in detail in Sec. II, represent length scales, /3i represent couplings and s encodes the time evolution 
of the characteristic mass scale of the model. 



defined in Eq. (38). The uncertainty in m(z), at any z-bin containing A'bin supernovae, is given by 

cTmiz)^ \ ^+dm^, (45) 

V ^^bin 

and we assume CTobs = 0.15. The systematic error, dm, is assumed to increase linearly with redshift: 

Am = bm , (46) 

^max 

(5m being the expected uncertainty and z,„ax the maximum redshift. In our analysis, we follow [90, 91] and assume 
(5m = 0.02 and z^^y^ = 1.7 for SNc from SNAP, plus low-z SNe from NSNF. Wc list the relevant parameters in Table 
II. The absolute magnitude, or the so-called nuisance parameter A^, is treated as an undetermined parameter in our 
analysis. 



IV. RESULTS 



Using the methods detailed in Section III, wc have evaluated the Fisher errors on the Modified Growth (MG) 
parameters {s, /3i, /?2, log(Af /Mpc^), log(A2/Mpc^)}, as well as the usual set of standard cosmological parameters. In 
Figs. 6-9 we plot the 68% C.L. contours for the MG parameters. Each of the ellipses is plotted after marginalizing 
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FIG. 7: The 68% confidence contours for the five MG parameters for Model II, based on f{R) with a different Compton scale 
(described in Sec. II) , as constrained by several different combinations of correlation functions from various experiments (shown 
in the legend). 



over all other parameters that were varied. The plots show contributions of different types of data: GC, GC+WL, 
and the combination of all data. We represent the contours from DES using different types of lines, while contours 
from LSST arc shown using different shades. 

Wc find that the five parameters are generally correlated and the degeneracy between log(A^/Mpc^) and 
log(A2/Mpc^) is strongly positive, while Pi and P2 have a strong negative correlation. The negative correlation 
between /3's is expected. One can enhance the growth by either raising (3i, which increases the effective Newton's 
constant, or lowering (32, which enhances the relative strength of the Newtonian potential ^' which drives the cluster- 
ing of matter (see Eq. (8)). This also explains the degeneracy between the two length scale parameters log(Af /Mpc^) 
and log(Ai/Mpc^). 

The full list of marginalized Fisher errors for the four fiducial models are listed in Table III for DES and Table IV for 
LSST. Comparing the constraints on the MG parameters from different datasets, we find that the constraining powers 
of WL, GC and WLxGC are comparable and that they provide much more stringent constraints than WLxCMB and 
GCxCMB. This is, in part, because there are many more different correlations between WL and GC than there are 
cross-correlations with CMB (see Fig. 2). Also, correlations of GC and WL with CMB suffer from a larger statistical 
uncertainty, since only the ISW part of the total CMB anisotropy is correlated with large scale structure. 

As expected, models with log(A2/Mpc^) = 4 (II, IV) are better constrained than models with log(A2/Mpc^) = 3 (I, 
III). In the former case, variations of the MG parameters give rise to larger changes in fi and 7 on linear scales. Also, 
the Pi parameters, which are directly related to the coupling a' in scalar-tensor theories, are much better constrained 
in Models III and IV, compared to Models I,II. This is because with s = 2, functions fi and 7 depart from their GR 
values of unity at much earlier times for the same scales (as can be seen in Fig. 1). Hence, with DES, and especially 
with LSST, one will be able to significantly reduce the volume of the allowed parameter space in Chameleon type 
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FIG. 8: The 68% confidence contours for tlie five MG parameters for Model III, based on a Cfiameleon model (described 
in Sec. II), as constrained by several different combinations of correlation functions from various experiments (shown in the 
legend) . 



models. Also, for all models we have considered, the Fisher errors are small enough to allow for meaningful consistency 
checks using Eq. (17). 

We can quantify the correlation between two randomly distributed parameters p, q by their correlation coefficient, 
which is defined as 

eb,.)^^^, (47) 

where Cov(p, q) and Cp, cr^ are the off-diagonal elements of the covariancc matrix and the variance of the parameters 
p,q, respectively. For example, with all data combined, ^[log(Ai/Mpc^), log(A2/Mpc^) ] ~ 0.9 and ^(/3i,/32) ~ —0.9 
for both models I and II. This strong degeneracy appears in both DES and LSST data, which reflects the degeneracy 
between ^ and 7. 

The scale-dependent growth pattern induced by modified gravity can, in principle, be partially degenerate with the 
bias factors of galaxies. Thus, in Fig. 10, we plot the correlation coefficients of the basic cosmological parameters 
(lower-left blocks), MG parameters (central blocks) and the bias factors (upper- right blocks) using model II and LSST 
as an example. In each panel, the off-diagonal blocks show the correlation among the basic, MG and bias parameters. 
In the case of GCxCMB, the MG parameters are correlated with bias for all the redshift bins, since GCxCMB is 
affected by modifications in our fiducial model over the entire redshift range of the LSST GC bins (see Fig. 3). Hence, 
the signal seen by GCxCMB can be either from modified gravity or from galaxy bias. For GC, the MG parameters 
are more strongly correlated with biases for redshift bins centered at z < 2 than for the higher redshift bins, because 
GC cannot detect the MG effects from our fiducial model at z > 2, as illustrated in Fig. 3. There is less degeneracy 
in the WLxGC panel than in the GC panel, since WL is independent of bias. Naively, one might expect much less 
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FIG. 9: The 68% confidence contours for tfie five MG parameters for Model IV, based on a Chameleon model with a different 
Compton scale (described in Sec. II), as constrained by several different combinations of correlation functions from various 
experiments (shown in the legend). 



degeneracy among bias and MG parameters after combining all of the data, including WL. However, since WL data 
contains degeneracies among basic parameters and MG parameters, the degeneracy cannot be effectively broken. This 
is what we see in the upper panel of Fig. 10. 

It is interesting to consider the case where we somehow have measured the bias independently. In Tables III and 
IV, we show the results with and without adding a strong prior on the bias parameters: a{bi) ~ 0.001. The results 
are listed in the A and B columns in Tables III and IV, and the correlation matrix is shown in the lower panel of 
Fig. 10. For the case considered in Fig. 10, the bias prior can help to improve the constraints on MG parameters, on 
average, by 19%, 42%, 30% and 31% for GCxCMB, GC, WLxGC and all data combined, respectively. Comparing 
these numbers with Fig. 10, one can see that the stronger the degeneracy between the MG parameters and the bias 
factors, the bigger improvement we can obtain by breaking the degeneracy with a bias prior. Moreover, thanks to 
this strong bias prior, the relic degeneracy among basic parameters, MG parameters and bias is effectively broken as 
shown in Fig. 10. 

We note the importance of not ignoring correlations between galaxy counts from neighboring bins. Including the 
galaxy cross-spectra is necessary for the estimation of covariances from the Fisher matrix (Eq. (31)), since bins do have 
an overlap. Its inclusion also helps to determine the bias and, hence, improve the constraints on the MG parameters. 
The cross-correlations between widely separated bins vanish and only the neighboring bins actually contribute. Still, 
without their inclusion, we find that the costraints on the MG parameters from LSST arc weakened by up to 40% in 
the case with no prior on the bias. Adding a strong prior on the bias nullifies the relative importance of the infomation 
contribution from galaxy cross-spectra. 

Even without the bias prior, we did not find significant correlations between the MG parameters and the basic cos- 
mological parameters, such as flbh^, ^ch^, etc., namely, the absolute value of the corresponding correlation coefficients 
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Model 


P fiducial 


1(7 errors 
GC+CMB+SN 
A B 


from Planck 
+WL 
A B 


, SNAP anc 

+WLxGC 
A B 


DES 

ALL 

A B 


I 


log(A?/Mpc^) 3 + log(3/4) 
Pi 4/3 

log(Ai/Mpc2) 3 
ft 1/2 
s 4 


3.8 1.9 
0.71 0.55 
3.8 2.1 
0.85 0.70 
6.1 2.6 


1.7 1.1 
0.46 0.34 

1.9 1.2 
0.59 0.45 

3.7 1.9 


0.99 0.71 
0.27 0.22 

1.0 0.75 
0.34 0.29 

2.1 1.4 


0.99 0.71 
0.27 0.22 

1.0 0.74 
0.34 0.29 

2.1 1.4 


II 


log(A?/Mpc2) 4 + log(3/4) 
Pi 4/3 

log(Ai/Mpc=^) 4 
132 1/2 
s 4 


1.7 0.58 
0.31 0.09 
2.1 1.2 
0.38 0.14 
3.5 0.95 


1.2 0.46 
0.22 0.08 

1.2 0.59 
0.29 0.11 

2.5 0.83 


0.68 0.35 
0.10 0.06 
0.70 0.41 
0.13 0.08 
1.2 0.67 


0.68 0.35 
0.10 0.06 
0.70 0.41 
0.13 0.08 
1.2 0.66 


III 


log(A?/Mpc^) 3 + log(8/9) 
Pi 9/8 

log(Ai/Mpc2) 3 
P2 7/9 
s 2 


2.3 0.83 
0.22 0.08 
3.5 2.5 
0.42 0.23 
4.2 1.2 


1.9 0.76 
0.16 0.06 
2.1 1.0 
0.31 0.14 
3.3 1.2 


1.2 0.62 
0.08 0.05 

1.3 0.78 
0.15 0.10 
1.8 0.95 


1.2 0.62 
0.08 0.05 

1.3 0.78 
0.15 0.10 
1.8 0.95 


IV 


log(A?/Mpc^) 4 + log(8/9) 
/3i 9/8 

log(Ai/Mpc2) 4 
ft 7/9 
s 2 


2.5 0.75 
0.10 0.03 
3.9 3.0 
0.22 0.11 
3.2 1.1 


1.4 0.66 
0.05 0.03 

1.6 0.99 
0.13 0.06 

1.7 0.89 


0.95 0.53 
0.03 0.02 

1.1 0.77 
0.08 0.04 

1.0 0.68 


0.94 0.53 
0.03 0.02 

1.1 0.76 
0.08 0.04 

1.0 0.68 



TABLE III: Error forecasts for parameters describing the modification of gravity based on the fiducial models considered in 
the paper, when DES is used for the GC and WL measurements. We list results for combinations of weak lensing (WL), 
galaxy number counts (GC), cross-correlations between lensing and counts(WLxGC), lensing and CMB (WLxCMB), counts 
and CMB (GCxCMB) and the combination of all of these (ALL). SNAP and Planck are always included in order to constrain 
the background parameters. We also illustrate the effect of adding a prior on bias parameters in columns A (no bias priors) 
and B (cr(bias)=0.001). 

are less than 0.3. 

Comparing the constraints from DES with LSST, we find that, typically, LSST can improve constraints on individual 
MG parameters by a factor of 3 or better. With all tomographic data combined (no priors used), the relative errors 
on MG parameters, which is the Icr marginalized error divided by the corresponding fiducial values, are summarized 
in Table V. There we can see that even DES will have significant power to constrain some of the MG parameters, 
giving 20% — 40% level constraints, while LSST can go below 10% levels. 

Given the constraints on MG parameters, one can reconstruct fj, and 7 using error propagation. The results are 
shown in Fig. 11 for Models I-IV. The inner shaded bands are the 68% C.L. regions from all combined data from 
LSST, and the outer lines show the constraints from DES. We see that even DES is able to put good constraints on 
/i,7 and test our fiducial models at high confidence level. In the plot, we see that the errors are generally smaller 
at small fc's and high redshifts. However, this does not really mean that we can constrain fi and 7 very precisely in 
that region. Instead, this is likely an artifact of our parametrization. One generally expects to find "sweet spots" , 
i.e. regions where the errors are small, close to the transition scale. This is analogous to the case of the dark energy 
equation-of-state w{z), where it is know that the number and the locations of "sweet spots" strongly depend on the 
parametrization of w(z) [92]. 

V. CONCLUSIONS 

In this paper, we have investigated the power of future tomographic surveys to constrain modifications in the growth 
of structure w.r.t. to that predicted in ACDM. Models of modified gravity, as well as models of coupled dark energy 
and dark matter, in general introduce a scale-dependence in the growth of structure and a time- and scale-dependent 
slip between the gravitational potentials. These modifications are expected to leave characteristic imprints on the 
observables, which in principle could be used to break the background degeneracy among different models of cosmic 
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Model 


P fiducial 


lo" erro 
GC+CMB+SN 
A B 


rs from Plane 
+WL 
A B 


k, SNAP and 

+WLxGC 
A B 


LSST 

ALL 
A B 


I 


log(A?/Mpc2) 3 + log(3/4) 
Pi 4/3 

log(Ai/Mpc2) 3 
P2 1/2 
s 4 


0.89 0.64 
0.27 0.20 
0.90 0.65 
0.33 0.25 
1.35 0.94 


0.39 0.27 
0.13 0.11 
0.39 0.28 
0.17 0.14 
0.73 0.54 


0.27 0.18 
0.08 0.07 
0.27 0.19 
0.10 0.09 
0.55 0.45 


0.27 0.18 
0.08 0.07 
0.27 0.19 
0.10 0.09 
0.55 0.45 


II 


log(A?/Mpc2) 4 + log(3/4) 

131 4/3 
log(Ai/Mpc2) 4 

132 1/2 
s 4 


0.40 0.21 
0.05 0.03 
0.42 0.24 
0.07 0.05 
0.57 0.32 


0.24 0.12 
0.03 0.02 
0.25 0.14 
0.04 0.03 
0.37 0.22 


0.15 0.09 
0.02 0.01 
0.15 0.09 
0.02 0.02 
0.24 0.18 


0.15 0.09 
0.02 0.01 
0.15 0.09 
0.02 0.02 
0.24 0.18 


III 


log(A?/Mpc2) 3 + log(8/9) 
Pi 9/8 

log(Ai/Mpc2) 3 
(32 7/9 
s 2 


0.58 0.29 
0.04 0.03 
0.70 0.42 
0.07 0.06 
0.68 0.31 


0.40 0.19 
0.02 0.02 
0.43 0.25 
0.05 0.04 
0.51 0.24 


0.25 0.14 
0.01 0.009 
0.27 0.16 
0.02 0.02 
0.32 0.20 


0.25 0.14 
0.01 0.009 
0.27 0.16 
0.02 0.02 
0.32 0.20 


IV 


log(A?/Mpc2) 4 + log(8/9) 
/3i 9/8 

log(Ai/Mpc2) 4 
/32 7/9 
s 2 


0.69 0.20 
0.03 0.007 
0.75 0.35 
0.06 0.02 
0.84 0.25 


0.38 0.17 
0.01 0.004 
0.39 0.21 
0.03 0.01 
0.40 0.18 


0.25 0.14 
0.007 0.003 
0.26 0.15 
0.01 0.007 
0.25 0.15 


0.25 0.14 
0.007 0.003 
0.26 0.15 
0.01 0.007 
0.25 0.15 



TABLE IV: Error forecast for the parameters describing the modification of gravity based on the fiducial models considered 
in Sec. II, when LSST is used for the GC and WL measurements. We list the results for combinations of weak lensing (WL), 
galaxy number counts (GC), cross-correlations between lensing and counts(WLxGC), lensing and CMB (WLxCMB), counts 
and CMB (GCxCMB) and the combination of all of these (ALL). SNAP and Planck are always included in order to constrain 
the background parameters. We also illustrate the effect of adding a prior on the bias parameters in columns A (no bias priors) 
and B (cr(bias)=0.001). 





Model I 


Model II 


Model III 


Model IV 


p 


DBS LSST 


DES LSST 


DES LSST 


DES LSST 


log(A?/Mpc2) 


34% 9% 


18% 4% 


41% 8% 


24% 6% 


Pi 


20% 6% 


9% 2% 


7% 0.9% 


3% 0.6% 


log(Ai/Mpc2) 


33% 9% 


18% 4% 


43% 9% 


28% 7% 




68% 20% 


26% 4% 


19% 3% 


10% 1% 


S 


53% 14% 


30% 6% 


90% 16% 


50% 13% 



TABLE V: Relative errors on the five modified growth parameters from all data combined. This table allows for a quick 
comparison between DES and LSST, and between different fiducial models. For more complete information see Tables III and 
IV. 

acceleration. It is therefore useful and important to explore to which extent the upcoming experiments will be able 
to detect and constrain modified growth patterns. 

We have used a five-parameter description for the rescaling of the Newton constant, ij,{a,k), and for the ratio 
of the metric potentials, ^{a,k), equivalent to that introduced in [32]. From the point of view of scalar-tensor 
theories (e.g./(i?) and Chameleon models) these parameters are related to the coupling in the dark sector and to 
the characteristic mass scale of the model. In these cases, the five parameters are not all independent; they need 
to satisfy the consistency conditions given by Eq. (17). This can be used to constrain the scalar-tensor models and 
potentially rule them out. In Sec. II we have described the /(i?) and Chameleon theories which we have used as 
fiducial models for our error forecasting. We have then studied in detail the constraints on the five parameters based 
on four fiducial models (two f{R) and two Chameleon) expected from Weak Lensing (WL), Galaxy Counts (GC), 
CMB and their cross-correlation spectra as seen by Planck, DES and LSST (additionally using Planck and SNAP to 
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FIG. 10: Correlation matrices for the basic, MG, and the bias parameters for the f{R) fiducial model with log(A2/Mpc^) = 4 
(Model 11). Parameters 1 to 5 are the basic cosmological parameters, namely, Qth^, 0.ch? , h, r and ris respectively; parameters 
6 to 10 are the MG parameters, namely, log(A?/Mpc^), log(A2/Mpc^), /32 and s respectively. Parameters 11 to 20 denote 
the 10 bias parameters for the LSST GC bins from low to high redshift. The positive and negative correlations are shown 
with shaded blocks with and without vertical lines. We compare the correlation matrices with(lower panel) and without(upper 
panel) imposing a bias prior of aibi) = 0.001. 



constrain the standard set of cosmological parameters). We have found that for scalar-tensor type models, DES can 
provide 20 — 40% level constraints on individual parameters, and that LSST can improve on that, constraining the 
MG parameters to better than 10% level. We have also found a strong degeneracy between the two functions /i and 
7, which was consistent with our expectations. Despite this degeneracy, however, the error ellipses are sufficiently 
tight for us to still find non-trivial bounds on the parameters of both functions. Overall, with DES, and especially 
with LSST, one will be able to significantly reduce the volume of the allowed parameter space in scalar-tensor type 
models. 

We have also found that the dilution of the constraints due to the degeneracy with linear bias is minimal when the 
full set of auto- and cross-correlation tomographic datasets is considered. 

The method we have employed relies on the choice of fiducial values for the MG parameters, and therefore is 
model-dependent. Nevertheless, the results we have obtained are encouraging as they show that upcoming and future 
surveys can place non-trivial bounds on modifications of the growth of structure even in the most conservative case, i.e. 
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FIG. 11: Reconstruction of the functions /i(a, fc) and 7(0, fc) from the constraints on the parameters {Af , s}, for Models I-IV. 
The inner (color-shaded) areas show the \a errors from LSST+Planck+SNAP, while the outer contours delimited by a solid 
line show the Ict errors from DES+Planck+SNAP. The line in the center denotes the fiducial model. 



considering only linear scales. This motivates us to apply the Principal Component Analysis (PCA) to the functions 
/X and 7. This method is more demanding computationally but, on the other hand, allows for a model-independent 
comparison of different experiments according to the relative gain in information about the functions. Furthermore, 
this method can point to the "sweet spots" in redshift and scale where data is most sensitive to variations in the 
functions /i and 7. This would be important information for designing future observing strategies. We present the 
results of the PCA analysis in a separate publication [35]. 

Another direction for future work is to study the effects of allowing for a small fraction of hot dark matter, such 
as neutrinos, and dynamical dark energy. One can, in principle, model the effects of dark energy perturbations and 
massive neutrinos by appropriately choosing the time-scale dependence of p. and 7. It is generally difficult to design 
a test that will definitively distinguish modified gravity from an exotic form of dark energy. However, it is clear that 
future data will have the potential to detect or, at the very least, significantly reduce the range of, departures from 
the ACDM model. 
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APPENDIX A: IMPLEMENTATION OF MODIFIED GROWTH IN CAMB AND SUPERHORIZON 

CONSISTENCY CHECK 

Following the notation in [74], we can map Eqs. (6) and (7) into synchronous gauge by using the following trans- 
formation, 

'^ = a + na, (Al) 
$ = 77 - Ha , (A2) 
a= {h + 6f])/2k^ . (A3) 

Now, in the synchronous gauge, Eqs. (6) and (7) become 

k\a + na) = ~^^lp^ , (A4) 

'n, = 7 ■ A5 

a + Ha 

Note that pA is gauge invariant. In the synchronous gauge, we have 

pA = p,„<5^ . (A6) 

and = Pb + Prnf^ni ~ P'o^h + Pc5% wlicrc subscripts m, b and c denote total matter, baryons, and cold dark 
matter, respectively, and the superscript s denotes the variables in synchronous gauge. Here and throughout, we 
assume that baryons comove with CDM in the late universe where we modify gravity. 

In addition, we have the energy- momentum conservation equations for cold dark matter and baryons [74] , 

(5^ = -i/i = 3?7- fc^a , (A7) 

-^'b = k , (A8) 

To solve the set of coupled differential equations (A4), (A5), (A7) and (A8) , we start by eliminating 77 in Eq. (A7). 
From Eq. (A5), we get, 

77 = 70; + (1 + ^)na . (A9) 

Taking the derivative w.r.t conformal time, 

7) = 7d + [(1 + -i)n + 7]d + [(1 + 7)7i + n^]a . (AlO) 
To obtain a and a we note that a is given by Eq. (A4), 

a = -H«-^MpA. (All) 

So, 

2 

a = -Ha -Ha- [pS^ + S^f, - Hp)] . (A12) 

Substituting Eq. (All) and (A12) into Eq. (AlO) and Eq. (AlO) into Eq. (A7) and doing some algebra, we can finally 
get the differential equation for 5*, which can be solved by using the a equation (All). In summary, after redshift 
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30, when gravity gets modified according to Eqs. (6) and (7), the equations we evolve are^ 



As _ Xs _ _ 
"c ~ "b ~ "m ~ 



-3aVm<[Ai7 + 7A - Wm(7 - 1)] + QM^ak^iH - H^) - 2M^k^a 



~n{k^a) 



~2Mf 



(A13) 
(AM) 



In our modified version of CAMB, we take the values of SI, (5,^ and a from GR at redshift 30 as the initial conditions 
to evolve Eqs. (A13) and (A14) at z < 30. 

Given a, d and d^-^, we have all the ingredients to calculate the observables, i.e the 's. For example, 



h = -2< , 

= 70; + aH{'j + 1) 
r, = (4 + k^a)/3 , 



a2/i 



+ * = {l + -f){d + na) , 

^J■ 

vj/ = + 7) + 7VI/ . 



' 2Af 2A;2 



(A15) 
(A16) 
(A17) 

(A18) 

(A19) 

(A20) 

(A21) 



As described in [42], if the theory is metric based, obeys causality in the infrared limit, and the energy-momentum 
conservation holds, then the evolution of super-horizon perturbations is uniquely determined if the relation between 
the two gravitational potentials is specified. Furthermore, in the long wavelength limit, fc — > 0, the ratio between 
the potentials will be scale-independent, i.e. 7 7(1)- It is important to test whether the equations (A13) and 
(A14) satisfy this consistency constraint on super-horizon scales. The evolution of the metric potentials on super- 
horizon scales is determined by the conservation equation for curvature perturbations [29, 42]. The relevant equations 
for the super-horizon dynamics are 



TTit 

H' 



IT 



(A22) 



where a prime, just in this section, denotes the derivative with respect to In a, and H = 

We can combine the equations (A22) into a second order differential equation for the potential 5* 



7 



El 
IF 



2_ 

7 











'IF J 


^1 



* = 0. 



(A23) 



We shall check whether equations (A13) and (A14) are consistent with Eq. (A23). In order to do so, we rewrite the 
equations (A13) and (A14) with independent variable In (a) and we combine them into a single first order differential 
equation for 



3ppa^j - 



2M2 fc2 



1 



ppa? 



* + aH'a - k^c 



(A24) 



In the super-horizon limit we can neglect the term 2Mpk^ ^ ippa?^ for any p and 7 which tend to a scale-independent 
finite function on large scales and obtain 



(A25) 



One should bear in mind that although cold dark matter and baryons evolve according to the same differential equation in the late 
universe, they have different initial conditions at redshift 30 when gravity gets modified. That's why we evolve them separately in 
MGCAMB. 
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Now, we can take the derivative of (A25) w.r.t. In a to obtain a second order differential equation for ^ 



7 



7" 7' 1 + 7' 



7 7 7 



aH" - aH' 



(A26) 



Finally, using the following equations 



a' H' ^ a 
aH'— = —--aH'- 

"f H ^ 7 

aH'a ~ -(1 + 7')* - 7*' 



(A27) 



it is easy to show that (A26) is equivalent to equation (A23). Therefore our equations give a consistent evolution on 
super-horizon scales. 



APPENDIX B: SCALAR-TENSOR THEORIES: MAPPING FROM THE EINSTEIN TO THE JORDAN 

FRAME 



The action for scalar-tensor theories in the Einstein frame reads 



S f 



d X\ —I 



(Bl) 



where Xi ^re the matter fields. (/) is a scalar field, and ai{(f>) represents the coupling of the scalar field to the i-th 
matter component. In what follows, we adopt the usual convention of indicating the Einstein frame quantities with 
a tilde. The inverse of the conformal transformation (11) maps the theory to the Jordan frame, where matter falls 
along the geodesies of the metric while the Einstein-Hilbert action is modified. Using this conformal map, we can 
easily determine the Jordan frame quantities in terms of the corresponding Einstein ones. At the background level, 
for the scale factor and the component of the energy-momentum tensor, we have 



e-2Ka.(0)^^ 



(B2) 



where U is the fluid 4-velocity and p and P are, respectively, the fluid energy-density and pressure. At the linear 
level, in the Newtonian gauge, we have 



= <If + 



(5 = (5 - 2Ka'^S4), 

V = V, 



Ka'fS(f) 



6P ^SP- 2Ka'/ 
a = a , 



(B3) 

(B4) 
(B5) 



where the prime denotes a derivative w.r.t. the field 



APPENDIX C: ANGULAR POWER SPECTRA 



Comparing predictions of cosmological models to observations involves projecting three-dimensional (3D) fields, 
such as the matter distribution or the gravitational potentials, onto two-dimensional (2D) fields measured by an 
observer looking at the sky. Consider, for example, a 3D field <%'(x, z), where x is the comoving coordinate and z is 
the redshift (which we also use as our time variable). It would be observed via its weighted projection - a 2D field 
X(n), which can be written as an integral along the line of sight: 

/•OO 

X(n)=/ dz Wx{z) X{hr{z),z) , (CI) 
Jo 

where n is a direction on the sky, r{z) is the comoving distance to a point at redshift z, and Wxiz) is a weight function, 
specific to the measurement, which selects a range along the radial coordinate that contributes to the observable X{h). 
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Ultimately, we are interested in the two-point correlation functions C^^ {9) = C"^^(|ni — n2|) = {X{hi)Y{n2)) 
between two fields - the auto-correlations (where X = Y) and cross-correlations (where X ^ Y). For convenience, 
we work with the Fourier transform of the sources 

Xirh,z)^ J ^Xik,z)e^-^\ (C2) 

The correlation function then becomes 

C^^(0) = J dz,Wx{zi) J dz2WY{z2) J ^e^'^ "^'-(^^) J l^e^k'.n.K^.) (^'(k, zi):y(k', z^)) . (C3) 

The temporal evolution of each Fourier mode depends only on the magnitude of its k-vector. Thus, we can write 
the source as a product of two factors. One contains the directional dependence at some early epoch z*, deep in 
the radiation era, when all modes of interest were well outside the horizon. The other factor represents the time 
dependence and hence controls the evolution of our source with redshift: 

X(k, f]) = z^)X(k, z) , (C4) 

where X{k,z^) = 1. Then, 

{X(k, zi)3^(k', Z2)) - X{k, zi)y{k', Z2) {X{k, z,)y(k', z,)) • (C5) 

At the high redshift z^,, when the modes are well outside the horizon and the modifications of gravity are negligible, 
the sources X(k, z*) and 3^(k, z^,) can be expressed in terms of the comoving curvature perturbation TZ [93]. Namely, 
one can always introduce coefficients cxn and cy-fi such that A'(k, z») = cxnT^O^, ^*) and similarly for y. Then, we 
can write 

{X{k, z,)y{k', z,)) = cxncyiz (7^(k, z,)7^(k', z,)) (C6) 

The homogeneity of space implies that 

(7^(k, z,)7^(k', z,)) = (27r)3,5(3) (k + k')Pn{k) , (C7) 

where Piz{k) is the primordial curvature power spectrum. One can also introduce the dimensionless spectrum, 
A|, = k^Pn/2-K'^, and write 

(A'(k,zOy(k',zO) =c;t7jcyK(27r)3<5(3)(k + k')^^^ (C8) 



The correlation function C-^^ {0) can then be written as 



C^y i^S) ^ cxncyn I dziWx{zi) l dz2iyj.(z2) / -^A|e''' ("i''i-"^'^^)<Y(fc, zi)3>(fc, zs) (C9) 



Integrating over the directional dependence of k, and defining R = \J f 'l + — 2rir2 cos 0, we obtain 

/■^* Z"^* t rik sin kR 

C''''{9)=cx^CY^ dz,Wx{zi)J^ dz2Wy{z2) J ^Al^^X{k,Zi)y{k,Z2). (CIO) 



We can further expand Eq. (CIO) in a Legendre series: 



c'^'^m^t.^-^crpdcose), (cii) 

where the first two terms in the series (the monopole and dipole contributions) are coordinate dependent and should 
vanish in the CMB frame of a homogeneous and isotropic universe. Expanding sui(kR)/kR in terms of spherical 
Bessel functions je, and comparing the result with Eq. (Cll), we can write 

Cr = ^^ J f A^/f(fc)/f(fc), (C12) 

where 

If{k) = cxn f ' dzWx{z)n[kr{z)]X{k,z). (CIS) 
Jo 

and similarly for ij . The expressions for /^'s for galaxy distributions in redshift bins, the maps of lensing shear for 
different bins, and the CMB temperature are given in Sec III A. One can then calculate various types of correlations 
between the different fields using Eq. (C12). 
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P flt,h'^ Qch^ h T Us 



Ap 0.005 0.005 0.03 0.02 0.005 



Model 


I 


II 


III 


IV 


P 




Ap 




log(A?/Mpc2) 


0.3 


0.4 


0.3 


0.4 


/3i 


0.1 


0.1 


0.1 


0.1 


log(Ai/Mpc2) 


0.3 


0.4 


0.3 


0.4 




0.05 


0.05 


0.1 


0.1 


S 


0.4 


0.4 


0.2 


0.2 



TABLE VI: The optimal stepsize for the cosmological parameters we use in this work. 



APPENDIX D: OPTIMAL STEPSIZE FOR FINITE DIFFERENCE DERIVATIVES 

For the calculation of the Fisher information matrix, we need to take numerical derivatives of our tomographic 
observables, say, Cf^ w.r.t cosmological parameters p. For small Ap, we can use the approximation 

dCf^/dp ^ [Cf^{p + Ap) - Cf^{p - Ap)]/2Ap (Dl) 

Note that Ap cannot be too big, which makes the right-hand-side of Eq. (Dl) deviate too much from the derivative, 
neither too small, which might give rise to numerical instability. To get a proper range for Ap is a bit of guesswork. 
One often has to do numerical tests until some range of Ap is found so that AC^'^ / Ap converges when Ap varies in 
this range. In Table VI we list the Ap we found and used in this work. 



APPENDIX E: RECONSTRUCTION OF /i AND 7 BY ERROR PROPAGATION 

Given the variances of the five parameters listed in Tables III and IV, namely, Var[log(Af /Mpc^)], Var[/?i], 
Var[log(A2/Mpc^)], Var[/32], Var[s], and the corresponding covariance matrices, one can estimate the variance of 
fi and 7 at given k and z by error propagation, fj. and 7 have the same functional form: 

/ - ii^, (El) 
■' 1 + t ' ^ ' 

t = Aik-Mpcfa", (E2) 

By perturbing Eq. (El) and (E2), we have 

A/ = [(/3-/)A< + iAA]/(l + <), (E3) 
At = i[lnlO- A(logA)+lna- As], (E4) 

Note that here we make A dimensionless for the ease of taking the logarithm later. So A is basically equivalent to 
log(A^/Mpc^) or log(A2/Mpc^). Taking the square of Eqs. (E3) and (E4) and plugging in the corresponding entries 
of the covariance matrix, we finally get 

Var[/] = (^)'{[lnlO • (/3 - /)]2Var[logA] + [Ina • (/3 - f)]^YaT[s] + Var[/?] 

+2{f3 - ff \nlO ■ Ina • Cov[logA, s] 
+2{(3 - /)lnlO • Cov[logA, f3] 

-h2(/?-/)lna-Cov[s,/3]|. (E5) 

where Gov [A, B] denotes the covariance of two random variables A and B, and the square root of Var[/] gives the 
errors on /x and 7. 
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